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ABSTRACT 

The  explicit  solutions  of  the  first-order  perturbation  equations 
for  a  one  dimensional  system  can  be  obtained  by  quadratures  provided 
that  the  zeroth- order  wave  function  does  not  have  any  nodes.  Methods 
of  Integrating  the  equations  are  developed  for  excited  states  where  the 
zeroth-order  wave  function  has  nodes.  One  of  these  methods  is  suitable 
for  obtaining  numerical  solutions* 
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THE  INTEGRATION  OF  THE  FIRST-ORDER  PERTURBED  WAVE  EQUATION 
FOR  EXCITED  STATES  OF  ONE  DIMENSIONAL  SYSTEMS 

Introduction 

Recently  there  has  been  a  renewed  Interst  in  perturbation  theory 

in  quantum  mechanics*  One  aspect  of  this  development^  due  largely  to 

Dalgarno  and  his  colleagues has  been  the  successful  calculation  of 

the  first-order  corrections  to  the  expectation  values  of  one-electron 

operators*  This  has  been  achieved  by  obtaining  the  exact  solutions 

of  the  first-order  wave  equations  rather  than  by  attempting  to  use  the 

very  awkward  representation  of  the  first-order  wave  function  in  terms 

of  the  spectrum  of  solutions  of  the  unperturbed  hamiltonian* 

2 

It  was  first  pointed  out  by  Young  and  March  that  the  first-order 
wave  equation  for  a  one  dimensional  system  can  always  be  integrated  by 
quadrature*  However  these  authors  only  discussed  the  ground  state  of 
the  system  in  which  the  wave  function  has  no  nodes.  In  the  case  of 
excited  states  singularities  occur  in  the  equation  owing  to  the  nodes 
of  the  unperturbed  wave  function*  The  object  of  this  paper  is  to  show 
how  the  singularities  may  be  handled  both  analytically  and  in  numerical 
work. 

First-Order  Wave  Equations 

For  a  system  with  the  hamiltonian 

H  =  H  +  V 
o 

the  perturbation  expansions  of  the  exact  wave  function  and  the 

exact  energy  E^  of  the  nth  state  may  be  written 


A*  Dalgarno^  "Quantum  Theory"  edited  by  D.  R*  Bates ^  (Academic  Press ^ 
New  York,  1961),  Vol.  I,  Chapter  5;  A*  Dalgarno,  "Advances  in  Physics" 
(Phil.  Mag.  Supplement)  II*  281  (1962);  C*  Schwartz,  Annals  of  Phys* 

2,  156,  170,  and  176  (1959). 

^  W.  H.  Young  and  N.  H.  March,  Phys.  Rev,  109*  1854  (1958). 
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♦n  ‘ 

Ao)  +*(i)  +  ... 

’n  ’^n  * 

(1) 

E  - 

E^°^  +  E^^^  +  E^^^  +  E^^^  +  •••  , 

(2) 

n 

n  n  n  n  * 

where  and 

^n  n 

are  the  nth  normalized  eigenfunction 

and  eigen- 

value  (assumed  non  •’degenerate)  of  the  unperturbed  hamlltonlan 


given  by 


(H  - 

o  n  ’n 


0  . 


The  equation  for  the  first-order  wave  function  is 


(H  -  +  (V  -  -  0 


(3) 


(4) 


and  therefore  the  first-order  energy  is  given  by 

,(1)  .  , 

■  •  <« 

The  second-  and  third-order  energies  are  given  in  terms  of  by 

,  (6) 

The  expectation  value  of  an  operator  L  in  the  nth  state  of  the 
perturbed  system  may  be  expanded  into  different  orders: 

n  n 


f 


(8) 


where 


3. 


(9) 


-  2  L 


(10) 


It  has  been  pointed  out  by  Dalgarno  that  in  cases  where  the  complexity 
of  the  perturbation  V  makes  the  first-order  equation  (4)  too  diffi¬ 
cult  to  solve  explicitly^  it  may  nevertheless  be  possible  to  solve  the 
first-order  equation 

(H  -  E^°^)X  +  (L  -  L  -  0  .  (11) 

o  n  '  ”  n  '  n'^n  '  * 


The  first-order  correction  (10)  may  then  be  written  in  the  form 


2 


(12) 


The  Integral  in  (12)  is  usually  tractable  even  if  V  contains  two- 
electron  terms  (such  as  that  the  correction  is  obtainable 

if  (11)  can  be  solved^  which  may  be  possible  if  L  is  a  one-electron 
operator.  The  mathematical  problem  is  therefore  to  solve  one-electron 
first-order  wave  equations  of  the  type  (11). 


Reduction  of  the  First-Order  Wave  Equation 

The  device  used  to  simplify  (11)  is  to  put 

■X,  ■  'A  •  <‘3) 

where  we  have  replaced  by  convenience,  and  to  note  that 

(H  -  *  -  (H  f  -  f  H  )6  ,  (14) 

o  n  n'n  on  n  o  'n 


so  that  (11)  may  be  written 
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where  the  square  brackets  denote  a  commutator.  This  means  that  if,  as 
usual j 


H  -  +  0  (16) 

o  o 

where  the  unperturbed  potential  energy  is  a  scalar^  then  (15) 

becomes 


If  this  equation  is  multiplied  by  ,  which  we  shall  assume  to  be 
real^  it  may  be  rearranged  to  give 

v.<*Jvy  -  ,  (18) 

which  is  equivalent  to  the  Wronskian  form 

-v*»  • 


Since  the  usual  boundary  conditions  are  that  and  vanish  at 

the  boundaries^  the  boundary  condition  on  f^  is  that 


0  on  boundaries  • 


(20) 


Solution  of  the  First*Order  Wave  Equation  for  One  Dimension 

For  a  one  dimensional  system  described  by  a  variable  x  confined 
to  the  interval  (a^b)  ,  equation  (18)  becomes 


(21) 


This  msy  be  integrated  directly  to  give 
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.2 

dx  “ 


M^(x) 


(22) 


where 


(23) 


and  the  constant  of  integration  vanishes  because  of  the  boundary 
condition  (20). 

The  same  treatment  may  be  applied  to  a  three  dimensional  system 
described  by  spherical  coordinates  in  which  the  operator  L  depends 
only  on  r  ^  and  the  unperturbed  wave  function  is  of  the  form 


(24) 


In  this  case  f^  depends  only  on  r  ^  and  the  equation  (18)  reduces  to 


4-(r^R^  ~)  -  2r^R  (L  -  L  )R 

dr'*  n  dr '  n'  n^  n 


(25) 


which  is  of  the  same  form  as  (21)  and  may  be  integrated  directly  from 
the  lower  limit  r  »  0  to  give 


2_2  . 


(26) 


where 


P^(r)  -  2  /'''s\(8)(L-L  )R  (8)<to  .  (27) 

n  wq  n  n  n 

(a)  Ground  state 

For  the  ground  state^  say  n  «  0  ^  in  which  the  unperturbed  wave 
function  ^^(3t)  does  not  possess  any  nodes  between  a  and  b  ,  (22) 
may  be  integrated  immediately  to  give 
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PX  M  (I) 

1  T2 - 

as  pointed  out  by  Young  and  March*  It  should  be  noticed  that  this 
Integration  Is  not  required  explicitly  If  *  to  be  used  In 

an  Integral  such  as  (12)*  This  Is  because 


-  if  f  <)  (V  -  V  )d)  dx 

O  J  0^0  o'^o 


may  be  transformed  by  partial  Integration  to  give 

t(1)  . 


where 


/b  df 

.  ir»< 


(x)dx  , 


■-/ 


b  M„(x)W  (x) 
o  o 


dx 


V*)  ■ 


In  particular  the  second^order  energy  becomes 


(29) 


(30) 


(31) 


(32) 


(33) 


The  only  point  requiring  special  discussion  In  the  radial  case 

when  R  (r)  has  no  nodes  Is  the  behaviour  of  £  (r)  at  the  lower 
n  n 

limit  r  •  0  a  If  &  (0)  #  0  .  then  f  (0)  Is  finite  as  long  as 

^  ^  -1  3 

LR  (r)  does  not  go  to  Infinity  faster  than  r  at  the  origin*  If 

^  j 

R  (r)  ss  at  the  origin,  then  f  (0)  Is  finite  as  long  as  LR  (r) 
n  R-1  ^  ® 

goes  to  zero  as  fast  as  r  at  the  origin.  The  formula  correspond¬ 
ing  to  (31)  In  this  case  Is 


-2 

If  LR^(0  goes  to  Infinity  even  as  fast  as  r  at  the  origin^  no 
serious  difficulties  are  encountered.  However^  In  this  case  fi|(0) 
might  be  Infinite*  For  example^  If  LR  (r)  goes  to  Infinity  as  r*^ 
at  the  origin^  then  ^^(0)  contains  a  Serm  k  ln(r)  * 
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where 


2  2 

rV(r) 


dr 


Qn(r) 


2  J  a\(a)Of-V^)tL^{s)da 


(3A) 


(35) 


(b)  Excited  states 


Consider  the  case  In  which  the  nth  unperturbed  eigenfunction 
*  ^n^^  ^  simple  zeros  at  a^^a2;** 

Then  If  the  first-order  wave  function  0( 


Is  well- 


(•,b). 

behaved,  the  function  have  simple  poles  at 

The  direct  Integration  of  (22)  to  yield  f  presents  difficulties 


•,a  In  the  Interval 
“  (1) 
n  "n 


because 


v> 


does  not  vanish  In  general  at  the  zeros 

(1) 


‘j 


Similarly 


the  Integral  for  L'  '  corresponding  to  (31) >  namely 


r(l) 


■I 


b  M^(x)W^(x) 


dx 


(36) 


Is  Infinite  If  vanishes  In  the  range  (a,b)  • 

To  avoid  the  difficulties  raised  by  the  singularities  in  (22)  and 
(36),  let  us  assume  that  <{l^,  f^  ,  and  are  analytic  functions  In 
the  vicinity  of  the  real  Interval  (a,b)  •  Then  equation  (22)  may  be 
written 


dz 


V^> 

♦n<*> 


(37) 


and  Integrated  from  a  to  x  along  any  contour  C  which  does  not 
pass  through  the  zeros  of  : 


aX  m  (z) 

f  (X)  -  f  (.)  +  *  -2 - dx  .  (38) 

“  "  J.  ♦’<«) 
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z-plane 


Figure  1, 


By  Cauchy's  theorem  C  may  be  any  of  the  contours  connecting  a  and 
X  shown  in  figure  since  the  Integrand  (38)  always  has 

double  poles  at  the  points  a^j***^a^  ,  and  these  do  not  contribute  any 
residues  to  a  contour  enclosing  them.  The  nature  of  the  poles  of 

follows  from  the  simple  character  of  the  zeros  of  and  the 

definition  of  by  equation  (23),  according  to  which  M;^(®j)  “  0  , 
and  therefore  In  the  vicinity  of  a^ 

M^(2)  -  +  O(z-aj)^  .  (39) 

The  aame  conalderatlona  apply  to  the  Integral  expreaalon  for  In 

n 

equation  (36),  and  the  correct  form  la 


M  (a)W„(r) 
n  n 


dz 


(40) 


where  the  path  of  Integration  C  from  a  to  b  avoids  the  double 
poles  at  • 

The  radial  case  can  be  dealt  with  In  the  same  way,  and  leads  to 
the  formula 


P  (z)Q  (z) 
n  n 

2  2 

z  R  (z) 
n 


dz 


(41) 
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It  follows  that  when  the  Integrations  in  (40)  or  (41)  ^  m  be  per¬ 
formed  analytically  It  is  unnecessary  to  consider  the  singularities 
explicitly^  since  formally  it  makes  no  difference  whether  one  integrates 
along  a  contour  with  complex  variable  z  or  along  the  real  axis  with 
variable  x  .  To  clarify  the  method  and  to  illustrate  the  power  of  the 
integral  formulae  for  excited  states^  it  is  useful  to  treat  a  simple 
example^  for  which  we  choose  the  polarizability  of  a  particle  in  a  one¬ 
dimensional  box* 


(c)  Example 

Consider  a  charged  particle  of  unit  mass  whose  position  is  described 
by  X  and  is  free  to  move  in  the  interval  (O^it)  •  The  eigenfunctions 
and  energies  are 

^_(x)  -  sln(nx)  ,  (n  -  .  (42) 

'u  f  1C  n 


How  suppose  a  perturbation  V  «  x  -  ic/2  is  applied^  corresponding  to 

the  energy  of  the  particle  in  an  electric  field.  The  polarizability 

(2)  (2) 

is  where  is  the  second-order  energy  of  the  nth  eigen¬ 

state  given  by  an  integral  similar  to  (40)^  namely 

f  j^W^(2)/(^^(z)j^dz  (43) 

^  o 

where  the  contour  G  from  0  to  ir  avoids  the  singularities  of  the 
integrand  at  kit/n  (k»l,2,*  •  •  ,n-l)  .  The  function  given  by 

w^(z)  -  -  (2-f)  j  .  (44) 

When  this  expression  is  substituted  into  (43)  and  the  integration 
performed  analytically  ignoring  the  singularities^  we  obtain  finally 


1  .15  2. 

— T  > 

24n  n 


(45) 


This  agrees  with  the  result  obtained  after  much  algebraic  labour  by 


10, 


summing  the  usual  spectral  formula 

CO 

,(2) 


E' 


I 


nk 


.1 

k»l  n  K 

(ife) 


Note  that  it  follows  from  (46)  that 

CO 


z 

n“l 


e<2)  -  0 

n 


which  can  be  used  to  check  equation  (45). 


(46) 


(47) 


Removal  of  Singularities 

When  the  integration  is  to  be  carried  out  numerically  by  a 
quadrature  designed  for  integration  along  the  real  axis^  it  is  necessary 
to  devise  ways  of  removing  the  singularities.  There  are  two  simple 
and  apparently  different  ways  of  doing  this^  which  are  however  related 
to  each  other. 


(a)  Removal  of  singularities  by  subtraction 


The 


function  f^  defined  by  (13)  has  n  simple  poles  at 


*  •  *a_ 


which  may  be  revealed  explicitly  by  defining  a  new  function 


J-1  J 


and  choosing  Che  so  that 


la  analytic  at  Che  points  a 


J 


(48) 


By 


aubstlCuClng  (48)  into  Che  differential  equation  (22)  we  get 


V> 


(49) 


The  coefficients 
singularities  of 


Cj^,*«*,c^  must  therefore  by  chosen  to  remove  the 

M  ;  that  is 
n  ^n  * 


II 


(50) 


Equation  (49)  may  now  be  Integrated  along  the  real  x-axia  to  give 

'xTH  (5)  -A  (I) 


8„(x)  -  g, 


nxfM  (J)  -A„(l)1 

[  Pld)  j 


dl 


(51) 


where 


n 

j-1 


(52) 


Similarly  equation  (40)  may  be  written  as  an  integral  along  the  real 
axis 

n 


s>  ■ 

•  -*  j-1  ^ 


(53) 


In  the  radial  case  equation  (48)  is  still  valid  with  the  coeffi¬ 
cients  given  by 


(54) 


The  equation  corresponding  to  (41)  then  becomes 

#f  r-w  -I-  -I  ^  ^-2* 

(1) 


where 


■  ^ 


dr  ,  (55) 


(56) 


(b)  Removal  of  singularities  by  lower  states 

An  alternative  method  is  to  make  use  of  the  n  eigenfunctions 
^o>***>^n-i  lower  states  of  the  unperturbed  system.  In  place 

of  (13)  we  write  the  first-order  wave  function  in  the  form 


n-1 


X„  -  8n^n  ■  Z  \K 

k-0 


12. 

(57) 


where  the  coefficients  b  ,  are  constants  to  be  chosen  so  that 

o'  '  n-1 

is  analytic  at  aj^,***,a^  .  By  substituting  (57)  for  into  (11) 

and  proceeding  as  before  we  get  in  place  of  (21) 

dg 

-  KK  ■  “  Z  ^<“1  -  V*A  •  <“> 

k»0 

which  on  Integration  gives 

dg  /»x 

* 

The  second  integral  may  be  evaluated  immediately  by  means  of  the 
unperturbed  SchrQdinger  equations  satisfied  by  and  ,  namely 


—r  +  2(=v  -  ^oA  "  ° 


dx 


A, 


(60) 


dx 


-  +  2(1I„  -  u  )6 


o  ^n 


By  multiplying  by  6  and  respectively  and  subtracting  we  get 

n  K 

d<b  d6. 

— (*  — -  -  6  -  2(E,  -  E  )  *,  6  .  (61) 

dx^\  dx  ^n  dx  ^  k  n'  “k^n  ^  ^  ' 


and  hence  by  integration 


♦kK-<'n<  ■  ■  V /' W. '■5  • 

a 

When  (62)  is  substituted  into  (59)  the  latter  becomes 

♦n  ^  -  Z  '■k'Vi  ■  • 

k-0 


(62) 


(63) 
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The  n  coefficients  have  now  to  be  chosen  so  that  the  two  terms 

on  the  right  hand  side  cancel  to  the  first-order  in  (x***j) 
points  aj  •  This  can  be  achieved  if  the  b^  satisfy  the  n  linear 
simultaneous  equations 
n-1 

Yj  ■■  •  (64) 

k-0 


It  follows  directly  from  (61)  that  the  derivative  of  the  right  hand 
side  vanishes  automatically  at  the  points  a^  ^  so  that  (63)  may  be 
integrated  along  the  real  axis  to  give 


where 


g^(x)  -  g^(a)  + 


n-1 


L  ^>^5)  i 


k-0 

The  corresponding  formula  for  then 


,(1) 


b,  V, 
k  kn 


where 


(65) 


(66) 


(67) 


(68) 


This  form  would  not  be  of  Interest  unless  It  was  possible  to  solve 

(64)  explicitly  for  the  bj^  In  a  large  class  of  cases.  The  class  Is 

that  In  which  the  unperturbed  functions  (^.  form  the  basis  for  a 

2 

Gauss- type  quadrature  formula  with  weight  factor  this 

case  quadrature  weights  A 


■j 


may  be  found  so  that 
n 


f  <)o(x)p(x)dx  -  ^  Aj^J(aj)p(aj) 

•  j-1 


(69) 


exactly  for  all  polynomials  p(x)  up  to  degree  2n-l  •  The  Aj  are 


given  by 


U 


4 

or  alternative  formulae  .  Applying  equation  (69)  to  p 
get 


(70) 


2 

o 


ve 


n 

^  Aj^^(aj)^j(aj)  -  (k,X-  0,l,-*-,n-l)  .  (71) 

j-1 

Therefore  by  multiplying  (64)  by  A  (^,  (a.)  and  summing  over  j  we  get 

J  K  J 

the  explicit  form 


b 


k 


j-1 


(72) 


The  class  of  functions  to  which  this  treatment  applies  Is  that  In 

which  Is  equal  to  times  a  polynomial  of  degree  k  In  x  (or^ 

more  generally^  in  a  transformed  variable).  Thus  It  Includes  the 
Hermlte  functions  for  the  simple  harmonic  oscillator  but  excludes  the 
Laguerre  function  solutions  of  the  hydrogen  atom. 

It  will  now  be  shown  that  In  the  cases  In  which  this  method  can 
be  applied  It  Is  equivalent  to  the  direct  subtraction  of  the  singular¬ 
ities.  Consider  the  expansion  of  the  function  <|)^(x) /(x-aj)  In  terms 
of  the  complete  set  of  unperturbed  functions  .  For  the  class  under 
consideration  only  the  n  lowest  functions  involved 

so  that 

n-1 

<^^(x)/(x-aj)  -  ^  Bj^(j)j^(x)  ,  (73) 

k-0 


where  the  B. 

k 


are  given  by 


/"’♦A 

x-a. 


dx 


(74) 


Z.  Ropal^  ''Numerical  Analysis'*  (Chapman  and  Hall^  London^  1955).^ 
chapter  7. 


Mow  the  quadrature  fomula  (69)  la  exact  In  thla  caae  ao  that 
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and  therefore 


n-1 

k-O 

(SO)  we  have 

n  n-1 

■  I  Z 

j-l  k-O 

n-1 

“  Z  • 

k«0 

where 

n 

■“k  ■  Z  ' 

j-l 

n 

■  Z  Vk'^j’^k'*)’'*"'*)’  ■ 

j-l 


^  . 

x-aj 

Hence  from  (48)  and 

4  8  -  X 

^n*n  ''  n 


(75) 


(76) 


(77) 


(78) 


by  equations  (40)  and  (72).  Thus  (77)  is  the  same  as  (57)  ^  which 
proves  that  in  the  restricted  class  of  cases  under  consideration,  the 
two  methods  are  identical;  the  same  arguments  apply  to  the  radial 
equation*  While  the  second  method  is  theoretically  interesting  in 


suggesting  that  quadrature  formulae  may  be  useful  in  perturbation 
theory^  the  first  method  is  quite  general  and  is  naturally  to  be 
preferred  in  numerical  applications* 
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